Color Confinement and Massive Gluons 



Masud Chaichian 1 and Kazuhiko Nishijima 1,2 

1 Department of Physics, High Energy Physics Division 
University of Helsinki 
and 

Helsinki Institute of Physics 
P.O. Box 9, FIN-00014 Helsinki, Finland 
2 Department of Physics, University of Tokyo, 7-3-1 Hongo, 
Bunkyo-ku, Tokyo 113-0033, Japan 



Abstract 

Color confinement is one of the central issues in QCD so that there are various 
interpretations of this feature. In this paper we have adopted the interpretation 
that colored particles are not subject to observation just because colored states are 
unphysical in the sense of Eq. (2.16). It is shown that there are two phases in QCD 
distinguished by different choices of the gauge parameter. In one phase, called the 
" confinement phase" , color confinement is realized and gluons turn out to be mas- 
sive. In the other phase, called the " dcconfinemcnt phase", color confinement is not 
realized, but the gluons remain massless. 



1 Introduction 

Strong interactions of hadrons are considered to be governed by QCD dealing with the 
interactions of quarks and gluons. It is a gauge theory based on the color SU(3) group 
and a strongly interacting particle belongs to an irreducible representation of this group, 
for instance, quarks and gluons belong to the color triplet and octet representations, 
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respectively. It so happens, however, that those particles that belong to non-singlet rep- 
resentations are not subject to direct observations and this property is referred to as color 
confinement. Thus both quarks and gluons are unobservable, and the only observable 
particles are color singlet composite particles called hadrons. 

Experimentally, we can study the properties of quarks and gluons only through hadrons. 
In order to study hadron-hadron interactions theoretically we shall first illustrate the prob- 
lem by a similar one in QED. The interaction between two charged particles is governed by 
Coulomb's law, but the interaction between two electrically neutral particles is represented 
by the van der Waals potential, 

V vdW {r)^r- & . (1.1) 

This shows that the electric fields generated by neutral systems penetrate into the vacuum 
without any sharp cut-off. 

In QCD we can elucidate the dynamical properties of color-neutral hadrons with ref- 
erence to dispersion relations. For the scattering of hadrons, for instance, they remain 
applicable provided that confinement excludes quarks and gluons from the physical inter- 
mediate states appearing in the unitarity conditions. This is precisely the condition for 
color confinement. As a typical example, let us consider nucleon-nucleon scattering, then 
the potential between them is given by the pole contribution in the crossed channels. The 
least massive hadron that can be exchanged between them is the pion, and the resulting 
interaction is represented by the Yukawa potential, 

Vy(r) oc (1.2) 

where /j, denotes the pion mass. 

In comparison with the van der Waals force in QED we recognize that the flux of the 
color gauge field emerging from color singlet nucleons cannot penetrate into the confining 
vacuum beyond a certain range thereby leaving no trace of long-range forces and that 
the penetration depth is given by the pion Compton wave length. Thus we recognize a 
similarity between the Yukawa mechanism for generating nuclear forces and the Meissner 
effect in the type II super conductors. This similarity strongly suggests that the vacuum 
allows penetration of the flux of the color gauge field generated by hadrons only by a finite 
length. 

Thus we many point out two salient features of QCD on the basis of the experimental 
properties of hadrons: i) color confinement, ii) finite penetration length of the fluxes of 
the color gauge field. 

The first one, color confinement, has been elucidated in [l]-[4] and we shall concentrate 
ourselves on the second problem in the present article. Actually, we interpret ii), finite 
penetration length, as an evidence for massive gluons and we shall show that realization 
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of color confinement as interpreted in [l]-[4] automatically leads to massive gluons. This 
is the main subject of this paper. 

In Sec. 2 we shall recapitulate the arguments for color confinement based on BRS 
invariance [7] and asymptotic freedom [5, 6]. In Sec. 3 we shall show that gluons turn 
out to be massive when color confinement is realized. In Sec. 4 we show that confinement 
is realized only in the gauges in which Z^ 1 = is realized. This is referred to as the 
confinement phase. On the other hand, gluons remain massless in the other gauges in 
which the above condition is not met, and this case is referred to as the deconfinement 
phase. 



2 BRS Invariance and Color Confinement 

In this section we shall recapitulate the essence of the interpretation of color confinement 
that has been developed in a series of articles [l]-[4]. 

In a covariant quantization of gauge fields introduction of indefinite metric is indis- 
pensable. Thus the resulting state vector space V involves unphysical states of indefinite 
metric and we have to find a criterion to select physical states out of V. For this pur- 
pose we employ the Lorentz condition in QED as a subsidiary condition, but it is more 
involved in non-Abelian gauge theories. In what follows we shall confine ourselves to 
QCD, and in order to fix the notation we start from its Lagrangian density in the metric 

= diag(l, -1, -1, -1) 

£ = C-inv + £gf + C-FP, (2-1) 

where 

C inv = ~ \FV ■ + WW Dp -m)fP, (2.2a) 
C gf = -A" ■ d^B + %B-B, (2.26) 

C FP = -i&*c- D^c (2.2c) 

in the customary notation [l]-[4]. We have suppressed the color and flavor indices in (2.2). 
The second Lagrangian density (2.2b) is the gauge-fixing term in which a denotes the 
gauge parameter and B the Nakanishi-Lautrup auxiliary field. The last one (2.2c) is the 
Faddeev-Popov ghost term, and the hermitian scalar fields c and c are anticommuting 
and are called Faddeev-Popov (FP) ghost fields. The local gauge invariance is respected 
only by the first term (2.2a) but not by the other two, (2.2b) and (2.2c), introduced for 
the purpose of quantization. The total Lagrangian is invariant, however, under the global 
BRS transformations [7] defined below. 

BRS transformations 
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Let us consider an infinitesimal gauge transformation of the gauge and quark fields 
and replace the infinitesimal gauge function either by c or c. They define two kinds of 
BRS transformations denoted by S and 5, respectively. 

5A^ = D^c, 5A /I = D^c, (2.3) 
Stf; = ig(c ■ T)ij), ~5ip = ig(c ■ T)tp. (2.4) 
where the matrix T is introduced in the covariant derivative of ip as 

= (fy - igT ■ A^, (2.5) 

For the auxiliary fields B, c and c local gauge transformations are not even defined, but 
their BRS transformations can be introduced by requiring the invariance of the local 
Lagrangian density, namely, 

SC = 5C = 0. (2.6) 

We shall not write them down explicitly, however, since they are not relevant to the 
following arguments. Noether's theorem states that the BRS invariance of the Lagrangian 
density amounts to two conserved BRS charges denoted by Qb and Q B . They satisfy 

H = i[Q B A] T , 5(f) = i [Qb.^] . (2-7) 

where we choose the — (+) sign, when the field 4> is of an even (odd) power in the ghost 
fields c and c. 

Equations of motion for the gauge field can be expressed with the help of BRS trans- 
formations as [l]-[4] 

d»F^ + gj u = i55A u , (2.8) 

where J v denotes the color current density and g the gauge coupling constant. It is worth 
noting that all three terms in Eq. (2.8) are divergenceless separately, in particular 

d u (i55A u ) = 0. (2.9) 
The BRS charges are hermitian and nilpotent, for example, 

Qb = Qb, Ql = 0. (2.10) 

The nilpotency implies introduction of indefinite metric and a physical state |/) is defined 
by the constraint 

Q B \f)=0, |/)€V. (2.11) 
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The set of physical states including the vacuum state |0) forms the physical subspace of 
V denoted by V phys , 

V ph ys = {\f):Q B \f)=0,\f)€V}. (2.12) 
Then the S matrix is BRS invariant and satisfies 

SS = i[Q B ,S] = 0, (2.13) 

so that the physical subspace V p h ys is an invariant subspace of the S matrix. 

Furthermore, we introduce a subspace of V called the daughter subspace Vd defined by 

Vd = {\f) : \f) = Q B \g),\g) eV}. (2.14) 

Then because of the nilpotency of Qb, Vd is a subspace of V p h ys , 

V d C V phys , (2.15) 
and we introduce the Hilbert space H by 

U = V phys /V d . (2.16) 

We may compare Vphys-, Vd and TL to closed forms, exact forms, and cohomology in Cartan's 
algebra and (2.16) may be called the BRS cohomology [3, 9, 10]. 

When the quark and gluon states are not physical, they are unobservable and hence 
confined. Thus the problem of color confinement reduces to demonstration of the condi- 
tions 

Q B \quark) ^ 0, Q B \gluon) ^ 0. (2.17) 

In a series of papers [l]-[4] it has been shown that the criterion for color confinement 
takes a simple form 

C = 0, (2.18) 

where the constant C is defined by 

d»(i5~8A a v {x),A){y)) = iS^CdjS^x - y) (j = 1,2,3). (2.19) 

Here and in what follows (• • •) denotes the vacuum expectation value of the time-ordered 
product. 

Furthermore, with the help of the renormalization group and asymptotic freedom [1]- 
[4], [10], it has been shown that the conditions (2.18) follows from 
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v = o, 



(2.20) 



where Z3 is the renormalization constant of the color gauge field. 

It has already been shown by Oehme and Zimmermann [11] without reference to per- 
turbation theory that the above condition is satisfied in the Landau gauge for Nf < 10 
on the basis of RGE and asymptotic freedom. It is true that the anomalous dimensions 
in RGE are evaluated in perturbation theory, but even then some results are valid be- 
yond perturbation theory thanks to asymptotic freedom. Depending on the nature of the 
objects RGE provides results at different levels of precision, namely, (1) approximate or 
semi-perturbative, or (2) exact. For instance, evaluation of the coefficient functions in 
the operator product expansion falls into the first category, whereas examination of some 
global properties of the theory falls into the second category. In fact, asymptotic freedom 
itself follows from the negative beta function evaluated in the lowest order, but this concept 
of asymptotic freedom is considered to be valid beyond perturbation theory. Likewise, the 
condition (2.20) for color confinement falls into the second category. This condition has 
been derived from unbroken non-abelian gauge symmetry and asymptotic freedom related 
to the high energy behavior of the quark-gluon system. When Z^ 1 is expressed as the 
integral of the Lehmann spectral function of the gluon propagator over the entire energy 
region, however, Eq (2.20) indicates a delicate balance of the contributions from all energy 
regions including both ultraviolet and infrared. In this sense this condition gives an in- 
direct constraint on the infrared behavior of this system. This should be contrasted with 
Wilson's area law in which the low-energy or long-distance behavior plays a dominant role 
in explaining quark confinement. It is also worth emphasizing that as a consequence of 
the above condition all the particles belonging to non-singlet representations of the color 
group are confined for the same cause. 

3 Pole Structure of Green's Functions and Massive Gluons 

In this section we shall show that gluons turn out to be massive when the condition for 
color confinement (2.18) is satisfied. For this purpose we start form Eq. (2.8). This 
equation has been given in the unrenormalized form so that we put the script (0) to 
unrenormalized expressions and rewrite (2.8) as 



Then we study its relationship to the renormalized version, and in order to facilitate 
understanding of the nature of the problem we start from its abelian version or QED, 
namely, 



(3.1) 




+ e Ji 



(3.2) 
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The multiplicative renormalization of fields and parameters relevant to this equation can 
be summarized as 

4°) = Z^A„ B^ = Z- 1 ' 2 B jM = J u , (3.3) 

e = Z 3 ~ 1/2 e, a = Z 3 a. (3.4) 
The renormalized version of Eq. (3.2) is given by 

d^F^ + eJ v = —d v B, (3.5) 

where 

eJ v = Z3- 1 [ e J„ + (1 - Z 3 )d l/ B] . (3.6) 

It so happens that both d^F^ u and d u B are multiplicatively renormalized but with different 
multiplicative factors. This mismatch forces to introduce operator mixing for renormal- 
ization [12]. Since both d^F^ u and d v B are finite operators so must be J u too. 

Essentially the same situation takes place in QCD, and the renormalized version of 
(3.1) is given by 

d^F- u + 9 J- = iSSA% (3.7) 

where J" is a linear combination of J" and i55A^. The space integral of Jg gives the color 
charge Q a . 

Q a = Jd 3 xJS(x), (3.8) 
satisfying the commutation relations of the color SU(3) algebra, 

Q a ,Q b ] =ifabcQ c . (3.9) 

With the help of Eq. (3.7) we can write down an equation for two-point Green's functions 
of the form 

(d x F^(x),A b u (y)) + (gJ-(x),A b u (y)) = {iSSA^A^y)). (3.10) 

In what follows we shall study the structure of the Fourier transforms of these Green's 
functions. 

Let F^ and G v be vector fields and introduce 

(F^(x),G„(y)) = j^pfrfke^-riT^ik), (3.11) 
and the Fourier transform of (F^, G u ) is denoted by 
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T^(k) = FT(F^G v ). 
Then T^ v can be expressed as a linear combination of two covariants: 

T,u(k) = -^T (fc 2 ) - - ^)T 1 (fe 2 ). 

Next we introduce two conditions: 
Condition 1): 

Assume d^F^ = and/or d^G^ = 0, then we get 

To(k 2 ) = To, constant. 

Condition 2): 

Assume in addition to the condition 1) that Q\0) = (0\Q = 0, where 



(3.12) 



(3.13) 



(3.14) 



Q = J d 3 xF (x), 



then we have 



(3.15) 



T n = 0. 



(3.16) 



Indeed, when this condition is not met we have a broken symmetry and the Nambu- 
Goldstone boson shows up in the form To ^ 0. 

All three terms in Eq. (3.10) satisfy the condition 1) and their Lehmann representations 
are given as follows. First, by taking account of the antisymmetry between subscripts A 
and u, we find 



FT(F X „ A v ) = -i{k x g» v - k,g Xu ) Lf^ + / dm 2 ^^ 



so that we obtain 



FT(d x F Xfl , A v ) = -R^- + (k 2 g, u " W / dm 2 ^g^. 
Then, thanks to the condition 2), we have in the absence of operator mixing 



FT(gJ^A u ) = -(feV " W/^ V-S+V 



Finally we have 



{k 2 g^ v - k^k v ) J rfm 2 i$i:- 



(3.17) 
(3.18a) 
(3.186) 

(3.18c) 



FT (iSSAfj,,A v ) = -C^- ~ \* y»u - j «»* W-m'+ie- 

In what follows we shall study the properties of these integral representations in more 
detail in QED and QCD. 

In QED we have R = C = 1, so that 
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FT(^F V ,^) = + (feV " W/^pzSe, ( 3 - 19 ) 

and 

FT(-^B,^) = -^, (3.20) 

so that we have 

FT{eJ„A u ) = (k 2 g, v - k,k v ) J dm 2 fc / ( ^V ( 3 - 21 ) 

The absence of the massless pole term in (3.21) reflects the fact that the conservation of 
charge is not broken. 

In QCD let us assume that the condition for color confinement (2.18) is satisfied, then 
we have 

FT(«5^, A v ) = -(fc V " W / dm 2 k .^l e , (3.22) 

which justifies the absence of the operator mixing in (3.18b). Since J M is a linear com- 
bination of and z<5<L4 M , we can form a linear combination of (3.18b) and (3.22) to 
find 

FT(gJ^A„) = -(k 2 g, u - k,k u ) J dm 2 J^^ , (3.23) 

where &2 is a linear combination of <T2 and 03. Then substituting (3.18a), (3.22) and (3.23) 
for Green's functions in Eq. (3.10) we find 

R = 0, (3.24) 

or 

FT(a A F v , A u ) = (k 2 g, u - k,k u ) J dm 2 1 ^±-, (3.25) 

Thus we find that there are no contributions to this two-point function from massless 
gluon states indicating the absence of the massless gluon. 

Then what would be the mechanism for generating the finite gluon mass? A possible 
candidate is the gluon condensate defined by 

{0\Afa)A b ft (x)\0)=5 ab K?0. (3.26) 

A BRS invariant modification of (3.26) has been considered by Kondo et alto introduce an 
effective gluon mass [13]. For phenomenological implications of such a gluon condensate 
as in (3.26), see e.g. [14]. In the Lagrangian density (2.2a) we find a term quartic in the 
gauge field and replacement of a pair of gauge fields by their vacuum expectation values 
(3.26) yields a mass term. By keeping only bilinear terms relevant to the description of 
the free massive gluon field, we find 
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£ 9 = \KA%, - \mlA»A^ (3.27) 

where 

A% = d,A a u - d v A«, (3.28) 

ml = glc 2 {G)K, (3.29) 

where C2(SU(N)) = N. This gluon condensate does not violate the color SU(3) symmetry 
nor the BRS invariance provided that the mass term in (3.27) is dynamically generated. 

4 Equivalence Class of Gauges 

In Sec. 2 we have introduced the QCD Lagrangian density that depends on the gauge 
parameter a. We shall study here how theoretical predictions depend on this parameter. 

In perturbation theory all the observable quantities are independent of the choice of 
the gauge parameter, but this is not the case in the non-perturbative approach as we shall 
see later in this section. Let us consider a class of Lagrangian densities {C a } representing 
a gauge theory such as QCD. Assume that all the members of this set are BRS invariant, 

SC a = 0, (4.1) 

and further that the difference between any two elements of this set is exact so that it can 
be expressed as the BRS transform of a certain operator A4, 

AC = C n - d = 8M, (4.2) 

then this set {C a } is called an equivalence class of gauges (ECG) in a loose sense. We shall 
show latter, however, that we have to introduce an additional condition for identifying the 
element of an ECG in the non-perturbative approach. 

Lagrangian densities corresponding to different choices of a in (2.2) belong to the same 
class in perturbation theory since we have 

A£ = ±(Aa)B-B = -i(Aa)5(c-B), (4.3) 

or 

M = -i(Aa)(c-B). (4.4) 

Now we introduce Green's functions in two gauges of the same ECG; then they are related 
to one another through the Gell-Mann-Low relation [15]: 
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(A( Xl )B(x 2 ) ■ ■ -) n = (A( Xl )B(x 2 ) ■ ■ ■ exp^iAS))!, (4.5) 
where A, B, ■■■ are local operators, and 

AS = J d A xAC = 5 J d 4 xM. (4.6) 
In particular, when all the local operators are BRS invariant, namely 

SA = SB = ■ ■ ■ = 0, (4.7) 

we obtain 

(A{ Xl )B{x 2 ) ■ ■ •)// = (A(x 1 )B(x 2 ) ■ • •)/, (4.8) 

by expanding the r.h.s of Eq. (4.5) in powers of AS. Since A, B,- ■ ■ are closed and AS is 
exact, we have 

{A,B,-..,(AST) = 0. (4.9) 

Since we are exploiting a series expansion in powers of AS, the proof of (4.8) is based on 
an assumed convergence of this series. 

In what follows we shall restrict ourselves to the Lagrangian density (2.2) and assume 
that all the quark masses have certain fixed values. Then a theory is characterized by a 
and g 2 . We shall plot (a,g 2 ) on a two-dimensional half-plane 

-oo < a < oo, g 2 > 0, (4.10) 

then each point on this half-plane corresponds to one theory. Now we introduce renormal- 
ization group (RG) with the generator of this RG given by [10] 

v = + Mi ~ 2a lv (g,a)^, (4.11) 

where fi denotes the renormalization point with dimension of mass and 7y is the anomalous 
dimension of the color gauge field. Since we have discussed applications of RG to Green's 
functions elsewhere [3, 10], we shall confine ourselves to the study of running parameters 
defined by 

g(p) = exp{pV) ■ g, 

a(p) = exp(pV) ■ a, (4.12) 
pZ(p) = exp(pV) ■ ii = e p p, 

where p denotes the parameter of RG. In what follows we shall concentrate our attention 
on g{p) and a(p). 
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When we increase p from to oo, the point(a(p),^(p)) moves along a line called the 
RG flow line (RGFL). Two theories corresponding to two points on the same RGFL are 
completely equivalent and physically identical as has been clarified in the applications of 
RG to Green's functions [3, 10]. The asymptotic limits of these running parameters are 
denoted by 

<?(oo) =5oo, a(oo) = a 00 , (4.13) 

and such a point represents a sink of the RGFL. Asymptotic freedom that has been 
assumed throughout this article is characterized by 

9oo = 0. (4.14) 

It should be emphasized that can assume only three alternative values depending on 
the choice of the starting point of the RGFL, namely 

«oo = -oo, 0, «o, (4-15) 
where ceo depends only on the number of quark flavors Nf in QCD 

oo = -i(13-|iV». (4.16) 
Without loss of generality we shall restrict ourselves to the following case: 

a > 0, or Nf < 10. (4.17) 
The value of depends only on the initial value of a, namely: 



Qtr 



Qo, for a > 0, 

0, for a = 0, (4.18) 
— oo, for a < 0. 



We recall the arguments on the ECG. Equivalence of two gauges characterized by 
(a±,g 2 ) and (a2,g|) was based on series expansion in powers of Aq, and when they 
belong to the same ECG so do also (o7i,gf) and (02,92)- The equivalence of the latter is 
based on the series expansion in powers of Aa(p). Now assume that a\ < and 02 > 0; 
however small the difference Aa = Q2 — ol\ might be, a\ tends to 00 and 02 to ao so that 
Aa(p) tends to 00. This fact casts a doubt on the convergence of the power series and we 
may conclude that these two points do not belong to the same ECG. Thus it is very likely 
that we have two sets of ECG defined on the two-dimensional parameter half-plane 

V(-oc) = {(a,g 2 ):a<0,g 2 > 0}, (4.19) 
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V(a ) = {(a,g 2 ):a>0,g 2 > 0}. 



(4.20) 



Both of them are two-dimensional domains, but in addition we have a one-dimensional 
line which forms the border between them, namely, 



L(0) = {(a ig 2 ):a = 0,g 2 > 0}. 



(4.21) 



Indeed, a(p) shows a discontinuity at a = when p is chosen sufficiently large as it is the 
case for a^. Then with the help of Eq. (4.9) of [4], i.e. 



G(pi;g,a,p) = exp 



Jo 



G(pi]g(p),a(p),p(p)) , 



where 7 denotes the anomalous dimension of the Green's function in question, we may 
conclude that Green's functions also develop a discontinuity at a = 0. This discontinuity 
is certainly reflected in the Green's functions, for instance, in the form of the residue of 
the massless pole of the gluon propagator studied in the preceding section. 

In this connection it is also important to recognize that particle masses defined as the 
pole positions of propagators are independent of p and thus are RG invariant as is clear 
from the above relationship. 

Thus it seems likely that we have three ECG specified by the value of a^. 

It has been shown already that the renormalization constant Z3 of the color gauge field 
is given by [3, 10] 



T-l 



a 



This constant is related to C in (2.19) through the formula [l]-[4] 

C{p) = a{p) - 2 f p °° dp'(j v (p) + j F p(p'))a(p 



x exp 



" 2 L P dp"j FP (p" 



(4.22) 



(4.23) 



where ~/fd is the anomalous dimension of the Faddeev Popov ghost fields and a(p) is given 
by 



a(p) 



(4.24) 



The validity of Eq. (4.23) has been tested in QED with the result that it reproduces C = 1 
correctly and the lack of asymptotic freedom in QED does not allow C to vanish. Hence 
charge confinement does not take place in QED. Therefore we have 



z, 1 



r 0, in D(-oo) and L(0), 
1 m V(a ). 



(4.25) 
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Combining (4.23) and (4.25) we find 



C = 0, in D(-oo) and L(0), 

C^O, inP(a ). K ' 

This amounts to the conclusion that gluons are massive in the confinement phase, T>{— oo) 
and L(0), whereas they remain massless in the deconfinement phase V(ao). 

The case of ao < is slightly more complicated, but essential features such as the 
existence of two phases are the same. 

5 Conclusions 

Theories described by the Lagrangian density (2.1) corresponding to different choices of 
the gauge parameter a belong to the same ECG in perturbation theory. This class is 
split into three, however, in non-perturbative approach depending on their convergence 
properties. Green's functions in QCD are expanded in double power series in g 2 and ag 2 , 
and the critical parameters governing the convergence of the expansion are the asymptotic 
values of g 2 , a, and ag 2 . The last one cannot be evaluated exactly, but is is still possible 
to judge whether it is equal to zero or non-zero. 

In order to clarify the differences among the three classes in the convergence properties 
we shall tabulate them in what follows: 

Case 1): q > or N f < 10. 



region 


{ag 2 )oo 


Z 3 


confinement 


gluon mass 


D(-oo): 







yes 


7^0 


D(0) or L(0): 








yes 


7^0 


D(a ) 





^ 0,oo 


no 






In this case D(0) or L(0) is a straight line a = and forms the border between D{— oo) 
and D(ao). 

Case 2): a < or 10 < N f < 16. 



region 


(ag 2 )oo 


Z 3 


confinement 


gluon mass 


D(-oo): 







yes 


^0 


D(a ) or L(q ): 





^ 0,oo 


no 





D(0) 





oo 


no 
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In this case D{olq) or L(a ) is a RGFL and forms the border between D(— oo) and 

Our naive belief in the independence of physics on the gauge parameter is no longer 
justified, but instead three phases show up. There are confinement and deconfinement 
phases, and in the former gluons turn out to be massive and in the latter they remain 
massless. 

Color confinement is one of the central issues in QCD so that there are various inter- 
pretations of this feature. In this paper we have adopted the interpretation that colored 
particles are not subject to observation just because colored states are unphysical in the 
sense of Eq. (2.16). Then the emergence of massive gluons is an inevitable consequence 
of this interpretation. 
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